Two-dimensional particle motion in a random potential under ac bias 



Dept. 



o 
o 

(N 

o 



S 

i 

C 

o 

(N 
> 

oo 

O 

o 

o 

-i— > 



c 

o 
o 



X 



Maxim A. Makeev/'B Imre Derenyi, 2 Q and Albert-Laszlo Barabasi'H 

' Dept. of Physics, University of Notre Dame, Notre Dame, IN 46566, USA 
of Biological Physics, Eotvos University, Pdzmdny P. stny. 1A, H-1 117 Budapest, Hungary 

(February 2, 2008) 



We study the Brownian motion of a single particle coupled to an external ac field in a two-dimensional random 
potential. We find that for small fields a large-scale vorticity pattern of the steady-state net currents emerges, 
a consequence of local symmetry breaking. In this regime the net currents are highly correlated, the spatial 
correlation function follows a logarithmic dependence, and the correlation length is of the order of the system 
size. For large external fields correlations disappear and only random net currents are observed. Numerical 
analysis indicates that the correlation length scales as a power law with both the size of the system and the 
amplitude of the ac field. 

PACS numbers: 79.20.Rf, 64.60.Ht, 68.35.Rh 



Driven diffusion in random media is a much studied prob- 
lem impacting on a number of various fields. By now, it is well 
established that the interplay between the annealed random- 
ness of the diffusion process and the quenched randomness of 
the media gives raise to unexpected scaling phenomena that 
have been extensively studied in recent decades HI SHI- Th e 
range of applicability of the problem of driven diffusive mo- 
tion in random media covers various fields, with relaxation 

rr 

phenomena in spin glasses |4], dislocation motion in disor- 
dered crystals [5], transport in porous media Jg], and turbulent 
diffusion 0] being but a few typical examples. In the presence 
of an external ac field, the motion of a particle is the result of 
the combined effect of thermally activated diffusivity and pe- 
riodic motion, induced by the coupling to the applied field. At 
time scales much larger than the period of the driving force, 
one would expect that the influence of the field is negligible 
and the dynamics of the system can be described by the clas- 
sical Brownian motion. Indeed, while the particle moves back 
and forth along the direction of the ac field, a stroboscopic 
view of the system, obtained by taking pictures only at times 
that are integer multiples of the external field period, would 
still show a randomly diffusing particle, almost as if the exter- 
nal field was absent in the system. In contrast to this intuitive 
picture, we demonstrate that the presence of an external ac 
field can fundamentally change the nature of the dynamics, 
when quenched randomness is present in the system. 

An indication, that an ac field bias may drastically influence 
the nature of diffusive motion in a random potential comes 
from the recent advances in the field of thermal ratchets deal- 
ing with one-dimensional (ID) systems, where the x —> —x 
symmetry is broken 0] • In equilibrium, a particle moving in a 
periodic asymmetric potential displays a simple diffusive be- 
havior. However, if the particle is also driven by an ac field, 
it drifts in the direction defined by the asymmetry of the po- 
tential, on time scales exceeding the period of the ac field. 
The appearance of this non-equilibrium steady-state net cur- 
rent is called the ratchet effect, and the average drift velocity 
of the particle is often referred to as the ratchet velocity |8]. It 
must be noted, however, that the ratchet velocity is expected 
to vanish when the particle is moving in a quenched random 



potential, since such a potential obeys inversion symmetry (in 
a statistical sense). In particular, the ratchet velocity in a sys- 
tem of identical energy wells separated by high enough energy 
barriers is always zero, irrespective of the distribution of the 
barrier heights. 

On the other hand, in two-dimensional (2D) systems any 
finite region of a random potential exhibits some degree of 
broken symmetry, and the particle can follow different tra- 
jectories during the two half periods of the ac driving. In this 
Letter, we show that the interplay between the local symmetry 
breaking in a random 2D potential landscape and the motion 
of the particle biased by an external ac field leads to the ap- 
pearance of a highly correlated steady-state net current field, 
which is characterized by a large-scale vorticity. We also in- 
vestigate the scaling properties of these fields. 

As a model system, we investigate the driven dynamics of 
a single particle moving in a random uncorrected potential 
on a 2D square Euclidean lattices (i, j) with varied number of 
identical sites N = L x L, where each site is connected to 
its nearest neighbors by bonds of unit length. Random poten- 
tial barriers are assigned to each bond and periodic boundary 
conditions are assumed. A schematic illustration of the sys- 




FIG. 1: Schematic illustrating the system under consideration. A 
particle on the site hops to its neighboring sites, with hopping 
rates, kj j, defined by both the potential barriers between these sites 
and by the ac field. 
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tern's geometry is shown in Fig. \l\ The particle, diffusing 
in the system, is also driven by an external ac field, applied 
along the x-direction. We define the local probability currents 
Jfj [Jfj] on site as the currents flowing over the po- 
tential barriers between the lattice sites and (i + 

+ 1)]. As shown in Fig. ^ a particle, located at site 
at an arbitrary moment of time t, can hop to any of the 
four nearest neighbor lattice sites by overcoming the potential 
barriers, E] p assigned to each bond. To simplify the nota- 
tions, we denote by ., E\ j, Ef, and Ef^ the local poten- 
tial barriers associated with the right E\ ^ = 

left E i,j = #(i,j)->(i-ij)> U P Efj = and 
down Efj = Efafi-tfaj-i) jumps, such that Ejj = E l i+l j 
and Efj = Efj +1 . The barrier heights are chosen to be 
random and quenched, with the spatial distribution given by 
— Eq + UrQj, where 7 stands for r, 1, u, and d; Eq is 

a constant; and rf- j = vl+i j an d Vtj = Vij+i are uncor- 
related random numbers uniformly distributed between and 
1 . Throughout the paper we measure the energies in units of 
IcbT, where /cb stands for the Boltzmann constant, and T is 
the absolute temperature. 

In addition to the thermally activated diffusion, the particle 
is also driven by an external ac force, 2<p(t), applied along the 
x-direction (see Fig.Q. Since the energy barriers are placed 
halfway between adjacent sites, the ac force modulates their 
hight by —(j>(t) in the positive x-direction and by +<f>(t) in 
the negative one. Thus, for large enough barriers, the hopping 
rates of the particle occupying site can be written as 

fc^.=^exp{-[t/77^. -0(t)]}, 

kl^veM-Pvlj + <!>(*)]}, 
ki,j= 1/e xp{- u '0i,j}, 

k$ ij =vexp{-Ur# ij }, (1) 

where v — vq exp{— Eq}, and the attempt frequency vq is 
also a constant. 

Next, we turn to an ensemble description, in which the 
probability of the particle occupying site is denoted by 
Pi.j, and the time evolution of this probability distribution is 
described by the following set of master equations: 

BP 

^ - - t x + j x - p> + jy (2) 

Q t - J i,j ^ J i-l,j ^ 

where 

jx yr p _ tJ p 

are the x and y components of the local probability currents 
Jij, as defined above. 

For simplicity, in the following we restrict ourselves to the 
consideration of symmetric square wave external fields, i.e, 
when the field <fi(t) alternates between +F and —F at con- 
stant time intervals Tp. We also assume that rp is much larger 
than the relaxation time of the entire system, estimated as 



Troiax ~ L 2 max(l / j) . In this case, for each half period 
of the ac driving the probability distribution Pjj and currents 
Jy relax to their steady state values: Py(+F), Ji,j(+F), 
and Pi j (-F), Jij (— F)), which can be determined from the 
stationary solution of the master equation for <f>(t) = +F 
and <f)(t) = — F, respectively. From direct analogy with the 
ratchet velocities, we can then define the net currents in the 
system as 

Ji,i = \[JiA+F) + JiA-F)\- (4) 

The appearance of non-zero net currents is a characteristics of 
systems with locally broken symmetry, and their magnitudes 
are given by higher than linear terms in the response functions 
{i.e., by F 2 and higher order terms) |9]. 

We solved the system of Eqs. Q-Q numerically by using 
the conjugate gradient method, modified for sparse matrices 
In Figs.|2l(a)-(d), we show the steady state local net cur- 
rent patterns obtained for a system with linear size L = 50, 
randomness parameter U = 0.5, and external field amplitudes 
F = 0.01, 0.05, 0.1, and 0.9. These figures offer a visual 
proof for the emergence of highly non-trivial steady-state net 
current fields, characterized by long-range correlations and 
large-scale vorticity. The vorticity structures are most appar- 
ent for small external field amplitudes (F U). As F in- 
creases, however, the correlations gradually disappear and the 
net current field converges to a nearly random structure. 

It must be emphasized that numerical solution of Eqs. Q- 
Q for arbitrarily small values of F and/or for arbitrarily large 
system sizes is not possible. This is due to the slow conver- 
gence of the conjugate gradient method in these limits and, 
most importantly, due to the growing impact of the random- 
ization, induced by the numerical errors |9]. 

To understand the nature of the observed correlations and 
investigate the scaling properties of the system, we compute 
the ensemble-averaged current-current correlation function, 
defined in the real space as 

(5) 

Here summation over all the lattice site indexes is im- 
plied, while summation over (i',f) indexes is performed only 
for lattice site pairs such that |rij — rv j/ \ = r, N r being the 
number of such pairs. The averaging < ... > was taken for 
various different realizations of the disorder {rj]j}. The cor- 
relation function defined this way is bounded from below by 
(in case of perfect correlation), from above by 4 (perfect anti- 
correlation), and takes the value of 2 for vanishing correlation. 

In Fig. |5J we show the normalized correlation function, 
C(r)/C(\), computed for different system sizes, and fixed 
external field amplitude, F = 0.01. As Fig. [3] demonstrates, 
the correlation function follows closely a logarithmic depen- 
dence on r, for small radial distances, and deviates from that 
as r grows beyond a certain value, with a saturation threshold 



3 



a) 



b) 



FIG. 2: Snapshots of the steady-state net current fields obtained for 
model systems with L — 50, U = 0.5, and different external field 
amplitudes: (a) F = 0.01; (b) F = 0.05; (c) F = 0.10; (d) F = 0.90. 



determined by the system size, L. This allows us to introduce 
a correlation length £ in the system, which we define as a char- 
acteristic length at which the correlation vanishes: C(r) = 2, 
that is, 1 - 0.5 C(r) ~< J(r hj )J(ri>j') > turns zero. 

On the basis of the above observations, we conclude that, in 
the most general form, the current-current correlation function 
is described well by the relation 



C(r) ~ log(r)/ 



(6) 



where f(x < 1) ~ const. Note that the correlation length 
in Eq. (|6} is dependent on both the amplitude of the exter- 
nal field and the system size; i.e., £ = £(F,L). Moreover, 
as our results demonstrate, the correlation length, £, scales 
with both of these quantities following power laws. Thus, for 
small amplitudes of the external field, the correlation length 
depends strongly on the system size, approximately following 
~ L a relation, where the exponent a = 0.983 ± 0.008 
(i.e., £ grows nearly linearly with L). On the other hand, for 
intermediate external field amplitudes, £ is found to be weakly 
dependent on L and to monotonically decrease with F, fol- 
lowing closely the inverse power law, £ ~ F~P, where the 
exponent is /3~ 0.78 ±0.01. The transition between the two 
regimes occurs at a certain field amplitude, F^(L), dependent 
on L. The above scaling laws can be collapsed into a single 
scaling relation: 



£~ F-Pg{FL a 'P), 



(7) 



where g(x) ~ x for F < F£(L) and g(x) ~ const for 
Ff(L) <F< F%(L). The significance of F§(L) will be dis- 
cussed later. In Fig. [4] we show the numerical data collapse, 
performed according to Eq. @. As one can observe, it pro- 
vides satisfactory results for small and intermediate FL a ^ . 
In general, three regimes of the scaling behavior can be dis- 
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FIG. 3: Normal-log plot of the current-current correlation function 
C(r)/C(l), computed for U = 0.5 and F = 0.01 for different 
system sizes. 
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FIG. 4: Scaling plot of the correlation length according to Eq. 
for different system sizes. 



current fields are observed to be largely independent of the 
particular realization of the disorder. The correlation length 
has been found to scale nearly linearly with the system size 
and display an inverse power law behavior with the external 
field amplitude. 

An intuitive explanation for the appearance of system-wide 
structures is that in 2D the particle has the opportunity to 
"choose" globally different paths when driven in the two op- 
posing directions. The superposition of these globally differ- 
ent flow fields results in large-scale flow patterns. This picture 
is also consistent with a slow (logarithmic) decay of the cor- 
relations and a correlation length being in the order of the sys- 
tem size. It is important to note that the above phenomenon 
cannot occur in ID, where all the energy barriers are in series, 
and no alternative trajectories are possible. 

Research at Notre Dame was supported by NSF-PHYS. 



criminated. For small F, the behavior is defined by the sys- 
tem size scaling and extends over the region of F values up 
to Fi(L). In the second regime, the scaling is dominated by 
the external field amplitude. This regime corresponds to the 
intermediate values of Ff (L) < F < F 2 C (L). Note that sub- 
stantial deviations from the scaling behavior predicted by Eq. 
are observed for large values of F {i.e., for F > jRj (£)). 
This behavior is an artifact of the discrete nature of the model 
system under consideration: since the particle is on the lattice, 
we cannot measure £ smaller than 1 . As a result, in the limit 
of large F, the correlation length does not follow Eq. Q down 
to arbitrary small distances, but rather saturates at £ ~ 1. 

In summary, we have studied the Brownian motion of a sin- 
gle particle on 2D Euclidean lattices of different sizes with 
quenched random potential, coupled to an external ac field. 
We have found that the interplay between the quenched ran- 
domness of the potential landscape and the external ac field 
bias leads to the emergence of large-scale vorticity patterns 
in the net steady-state currents. These currents are found to 
be strongly correlated, with the two-point correlation function 
following a logarithmic dependence. The properties of the 
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